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Electroweak phase transition by four dimensional simulations* 
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The finite temperature phase transition in the SU(2)-Higgs model at a Higgs boson mass Mb — 35 GeV is 
studied in numerical simulations on four dimensional lattices with time-like extensions up to Lt = 5. T c /Mh 
is extrapolated to the continuum limit and a comparison with the perturbative prediction is made. A one-loop 
calculation to the coupling anisotropics of the SU(2)-Higgs model on lattices with asymmetric lattice spacings is 
presented. Our numerical simulations show that the above perturbative result is applicable in the phenomenolog- 
ically interesting parameter region. 



1. Introduction 

At high temperatures the electroweak symme- 
try is restored. Since the baryon violating pro- 
cesses are unsuppressed at high temperatures, 
the observed baryon asymmetry of the universe 
has finally been determined at the electroweak 
phase transition. Due to the bad infrared fea- 
tures of the theory, the perturbative approach 
predicts 0(100%) corrections and breaks down 
in the physically interesting region, Mji > 65 
GeV . In recent years lattice Monte Carlo sim- 
ulations have been used in three (e.g. and 
four dimensions (e.g. ||) to clarify the details of 
the phase transition. Since the bad infrared be- 
haviour is connected with the bosonic sector of 
the Standard Model, in both cases the fermions 
are omitted and the SU(2)-Higgs model is stud- 
ied. The comparison of the two approaches gives 
not only an unambigous description of the finite 
temperature electroweak phase transition, but an 
understanding of the non-perturbative features of 
the reduction step, too. 

In this talk two selected new results of the four 
dimensional formulation are discussed. In Sec- 
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tion 2 the SU(2)-Higgs model at a Higgs boson 
mass Mh — 35 GeV is studied. Our four dimen- 
sional, finite temperature lattices have time-like 
extensions up to L t — 5. We extrapolate T c /Mh 
to the continuum limit and compare the result 
with the perturbative prediction. In Section 3 we 
present the one-loop calculation to the coupling 
anisotropics of the SU(2)-Higgs model on lattices 
with asymmetric lattice spacings. We test the 
results by numerical simulations and show that 
the perturbative result is applicable in the phe- 
nomenologically interesting parameter region. 

The simulations have been performed on the 
APE (Alenia Quadrics) computers at DESY-IFH 
and on the CRAY-YMP at HLRZ Julich. 



2. Critical temperature for Mjj ~ 35 GeV 

The details of the results presented in this Sec- 
tion can be found in 

Our studies at Higgs boson masses ~ 20, 35, 
50 GeV showed that the perturbative predictions 
are in quite good agreement with the lattice re- 
sults for the jump of the order parameter, latent 
heat and interface tension jl]||,||]. Unfortunately, 
the non-perturbative results are obtained on lat- 
tices with temporal extensions L t = 2, 3 and they 
have typically O(10%) errors, thus the compar- 
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ison with continuum perturbative predictions is 
quite difficult. 

The situation is much better for the critical 
temperature. Its value is 1/L t in lattice units. 
In order to extract the dimensionless T c /Mh one 
needs the transition point at a given Lt and the 
T = mass of the Higgs boson at that point. 
Both of them can be measured quite precisely. 
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Figure 1. Numerical results for T c /Mh versus 
(aT c ) 2 = L^ 2 . The straight line is the extrapolation 
to the continuum value shown by the filled symbol. 
The dashed lines are the perturbative predictions at 
order g 3 (upper) and g 4 (lower), respectively. Rhw 
gives the mass ratio of the Higgs and W boson masses 
and g\ is the renormalized gauge coupling. 



Fig. 1 shows our results. The non-perturbative 
masses are determined by a careful extrapolation 
to infinite volumes. As it can be seen, the T c /Mjj 
predictions of the one-loop and two-loop pertur- 
bation theory almost coincide. The errors of the 
lattice data are dominated by the uncertainties 
in the critical hopping parameter. The extrap- 
olated continuum value is T C /M H = 2.147(40). 
The value at L t — 2 is about 5% smaller. This 
relatively small deviation is better than the ex- 
pectation based on lattice perturbation theory. 
The value of T c /Mh extrapolated to the contin- 



uum limit differs by about three standard devia- 
tions from the two-loop perturbative result. This 
is under the assumption that L t — 2 can be in- 
cluded in the extrapolation, which is supported 
by the good quality of the fit (% 2 ~ 1). 

3. The SU(2)-Higgs model on asymmetric 
lattices 

The details of the results presented in this Sec- 
tion can be found in 

For larger M H (e.g. M H = 80 GcV) the elec- 
troweak phase transition gets weaker, the lowest 
excitations have masses small compared to the 
temperature, T . From this feature one expects 
that a finite temperature simulation on isotropic 
lattice would need several hundred lattice points 
in the spatial directions even for L t = 2 tempo- 
ral extension. This difficulty can be solved by 
using asymmetric lattices, i.e. lattices with dif- 
ferent spacings in temporal (at) and spatial (a s ) 
directions. The asymmetry of the lattice spac- 
ings is characterized by the asymmetry factor 
£ = a s /a t . The different lattice spacings can be 
ensured by different coupling strengths, (f3 s and 
flt for space-space and space-time plaquettes; k s 
and n t for space-like and time- like hopping terms) 
in the action for different directions (c.f. ||]). The 
anisotropics jp = fit/ Ps and 7^ = Kt/ K s are func- 
tions of the asymmetry £. On the tree-level the 
coupling anisotropies are equal to the lattice spac- 
ing asymmetry; however, they receive quantum 
corrections in higher orders of the loop-expansion. 
On the one-loop level one gets 
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where g and A are the bare gauge and scalar cou- 
pling, respectively. 

In general, the determination of (£) and 
7 K (£) should be done non-perturbatively. This 
can be achieved by tuning the coupling strengths 
requiring that the Higgs- and W-boson correla- 
tion lengths in physical units are the same in the 
different directions. 

We have used the above idea in perturbation 
theory The one-loop results are summarized on 
Fig. 2. As it can be seen only cp(^) and c K (£) 
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are given. The functions fya(£) and 6 K (£) vanish 
in this order. Our calculation is the extension of 
to a gauge-Higgs model. Omitting graphs con- 
nected with the scalar sector, the result of can 
be reproduced (the function cp{t) of the present 
paper corresponds to c T (£) — c CT (£) of ref. 0). 

It is necessary to check the validity of the 
perturbative result by non-perturbative methods. 
M H ~ 80 GeV has been chosen, thus M H ~ 
Myy at zero temperature with typical correla- 
tion lengths of 2-4 in lattice units. We have 
performed simulations with two different sets of 
anisotropies. First 7/3 = j K = 4 (tree-level). 
Then 7/3 = 3.8, 7^ = 4. As it can be seen on 
Fig. 3 in the first case £w > £h (£« = Mi(space — 
direction) /Mi(time — direction), where i denotes 
the W or Higgs chanel). In the second case 
£,w < £,h- Rotational invariance is restored if 
£,w = £,h ■ Our linear interpolation is shown by 
the two solid lines. The "matching point" is in 
complete agreement with the perturbative predic- 
tion (full triangle) . In addition we have simulated 
at the perturbatively predicted couplings. At this 
third point the rotational invariance is restored, 
thus £w = £,h within errorbars. 

The non-perturbative features of the theory ap- 
pear at high temperatures. Since the theory is de- 
fined at T = 0, where it is weakly interacting, the 
good agreement between perturbation theory and 



non-perturbative results is not surprising. The 
perturbative predictions open the possibility to 
study the phase transition for higher Higgs boson 
masses and to vary the space-like and time-like 
lattice spacings separately. 

Two of us (F. Cs. and Z. F.) were partially sup- 
ported by Hungarian Science Foundation grant 
under Contract No. OTKA-F1041/3-T016248/7. 
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